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This invention relates to guidance systems for 
radar homing missiles and has an important application 
to anti-aircraft missiles employing a guidance system 
which is capable of acting in an active mode. An 
active mode is a mode in which the target is illuminated 
by radiation from the missile and in which reflected 
radiation from the target (termed "skin echo") is 
received by the missile during its flight, the missile 
having been launched on a course to intercept the target 
by a guidance system attached to a launching aircraft 
or platform* The reflected radiation is processed to 
obtain the bearings, azimuth and/or elevation, and the 
speed or range of the target so that the missile can 
follow the changes of direction and speed of the target • 

An object of the present invention is the 
improvement of terminal accuracy in a pulse doppler 
radar seeker. 

According to the present invention a missile 
guidance system Qomprises means for transmitting pulsed 
radiation to illuminate a target, a multi-element 
aerial arrangement arranged to respond to reflected 
radiation and provide a plurality of correspondingly 
pulsed output signals, receiver means for deriving 
from the aerial output signals sum and difference signals, 
together providing an indication of the presence and 
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position of a target, the receiver means including means 
for tracking a target in doppler frequency or in angle 
and adjusting the parameters according to the conditions. 

The relevant features of one embodiment of a 
pulse doppler radar seeker will now be described , by- 
way of example, with reference to the accompanying 
drawings, of which :- 

Figure 1 is a basic block diagram; 

Figure 2 is a block diagram of a doppler 
tracking loop applied to the arrangement of Figure 1; 

Figure 3 is similarly a block diagram of an 
angle tracking loop applied to the arrangement of Figure 1, 

Figure 4 is a table illustrating a combination 
of target confirmation tests, 

and Figure 5 is a block diagram of an angle 
tracking and servo-control system. 

In the guidance system shown in Figure 1 there 

are three front mixers 9, 10 and 61 to which the sum 

signal S, the azimuth difference signal D and the 

a 

elevation difference signal D are applied. The three 
mixers 9, 10 and 61 feed range gates 23 in the three 
channels 19, 20 and 21. A further channel 62 to which 
the sum signal S is fed by way of mixer 9 has a split 
range gate 63 which is initially opened for only the 
first half of a range gate period and is inverted in 
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phase (i.e. given a gain of -1 to provide an output signal 
at an amplifier 64. The output signal from the amplifier 
64 is applied to a Fast-Fourier- Transform filter bank 
(FFT) and thence into a range gate control circuit 35 
which adjusts the opening times of all the range gates 
23 to a time delay T as required to track the target in 
range. 

At the output of the head amplifiers 13, 14 

and 15 in the channels 20, 19 and 21, the difference 

signal D is added to a fraction k of the sum signal 

S (produced by a voltage controlled attenuator 66) to 

produce an output signal which is fed to the I.F. 

amplifier 25. The output signal from the I.F. amplifier 

25 is fed to an electronic angle tracking phase-sensitive 

detector, to be described later, the output of which is 

zero when the difference channel 20 has zero signal. 

This occurs when k e S = D e » If the output of the 

difference channel 20 is not zero it passes the signal 

to an integrator 'which controls the value of k until 

e 

the relationship k S * D is satisfied. Since the 

target angle off aerial boresight is given by D /S then 

e 

k e equals this angle. 

Similarly at the output of the head amplifier 
14 the difference signal D fl is added to the sum signal S 
attenuated by an amount k a in a voltage controlled 
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attenuator 67 to produce an output signal which is fed 

to the I.F. amplifier 25. The output signal from the 

I.F. amplifier 25 is fed to an electronic angle tracking 

phase-sensitive detector, the output of which is zero 

when the difference-channel 19 has zero signal' in it. 

If the output of the difference channel 19 is not zero 

it passes the signal to an integrator which controls the 

value of k Q until the relationship k S = D_ is satisfied. 

Since the target angle of aerial boresight is given by 

Da/S then k Q equals this, 
a 

The Fast-Fourier-Transf orm circuits provide 
outputs giving the complex amplitude in each of the FFT 
filter •bins 1 (typically 16, 32, 64-, 128 or 256 •bins 1 
according to the bin width control setting.) The filter 
bin width is equal to the total FFT bandwidth divided by 
the number of filter bins and the output data rate is 
proportional to the bin width. Figure 2 shows how the 
sum FFT is used to control the doppler tracking loop 
and for target detection and conf irmation. 

Target detection is carried out on. all the 
power outputs (i.e. the square of the sum of the real 
and imaginary components) of the FFT with the exception 
of the first and last N L where N L = 9, 5, 3, 2, 2 for 
256, 128, 67, 32, 16 bin FFTs respectively. The power 
output is divided by the average power of all the bins 
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excluding the first and last N L , and then this value 
is compared to a threshold (T D ) which is determined so 
as to allow a certain number of crossings if the input 
is pure thermal noise (false alarms). The technique 
used is to compare each bin, scaled by the average 
power, with the maximum of the previous bin and a 
threshold T^. The maximum bin at the end of the process 
which also exceeds T D is taken as a target alarm. If 
no bin exceeds the threshold then no alarm is found 
and the detection process is applied to the next set of 
data. This detection method would be modified in 
practice to allow multiple alarms. 

If a target alarm is found its bin number is 
output so that the doppler tracking loop and confirmation 
processes can act on the correct frequencies. An 
additional calculation is performed in the detection 
process. The contents of the alarm bin are reduced by 
the threshold. This is done because the probability 
distribution of the maximum bin which exceeds a given 
threshold is approximately a constant (T^) plus a 
Rayleigh distribution (provided T D is > log e N). By 
subtracting the threshold value, the initial alarm 
can be treated in the same way as subsequent power 
outputs from the bin and this simplifies the confirmation 
system. 
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The calculations performed for each bin power 
p. are thus 




i-Pi/ > *3 (1) 



If r ± > T D and also > r c max then set r c max = r ± and 
store i as the alarm bin number. r c max is the maximum 
value of 'so far 1 . The stored value of i which 
remains when all bins have been considered is then the 
target alarm bin. This initial setting, or estimate, 
of the target frequency in Figure 2) is then employed 
to close the doppler tracking loop and initiate the 
target confirmation process. 

In order to get a high probability of target 
acquisition with a low probability of false confirm 
it is necessary to incoherently sum over a large number 
of FFTs. To make this possible target tracking is 
necessary, since 'either the target must be kept in a 
fixed bin or the alarm bin must be known in some other 
way. For this reason a two stage acquisition process 
is used. The first stage which has already been 
described detects the whereabouts of a likely target, 
the second, confirmation, ensures that this detected 
signal is indeed a target with a high degree of certainty, 
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The method of operation of the confirmation 
process is as follows:- the 'alarm' bin is designated 
by the doppler tracking algorithm, except immediately 
after the first detect when the detection subroutine 
provides this information. As in the detection algorithm, 
the equivalent of Equation (1) above is used to provide 
a scaled bin output. (If the alarm bin is out of the 
range N L +1 to N-N L the scaled bin is set to zero). Then 
this scaled 'target bin' on one FFT is added to the same 
quantity on the successive sample. This summing continues 
until a certain number of FFTs ('NFF') have been added 
together. At this point the sum is compared to two 
thresholds, an upper and lower threshold. If the sum of 
the bins is greater than the upper threshold then a 
•confirm' flag is set, if the sum is less than the lower 
threshold the suspected target is rejected. If the sum 
is in between then the flag remains in a 'don't know' 
state. If the target is not rejected a further set of 
'target' bins from 'NFF' FFTs are added to the first set 
of 'NFF' and the sum of 2' NFF * FFTs compared to two new 
thresholds, the upper one again denoting 'confirm' and 
the lower one 'reject'. After a certain number of 
stages ('NST') (provided an input to the program) the 
upper and lower thresholds are made to coincide so that 
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a definite decision is made. This is called the 
increasing sum test. 

At this point a moving average of NSTxNFF FFTs 
begins and the criterion for confirmation is that the 
moving average should be above a threshold (normally 
taken to be the last stage threshold) and that in 
addition the sub-sets of •NFS 11 FFTs should not reject. 
This latter test is put in to ensure reasonably rapid 
response if the target should suddenly disappear for 
some reason. Two tests are therefore performed. The 
first is to see if the moving average is above a simple 
threshold or not, and the second is the test on the 
last NTT (or 2NFF etc) FFTs compared to two thresholds. 
The combination of the two tests lead to one output and 
these are shown in the table of Figure 4. 

In addition a certain amount of re-initialisation 
takes place. If the overall confirmation goes off 
•confirm* the moving average is stopped. Also if the 
increasing sum test rejects or confirms the increasing 
sum stage is re-initialised (otherwise it continues 
by adding in the next , NFF' FFTs), and each * NFF 1 FFTs 
are tested separately. 

When target confirmation is complete the 
'confirm 1 output is used to enable the angle tracking 
loop and the missile guidance loop. 




The Doppler Tracking system employed is to 
allow doppler tracking to take place both within the 
digital processor of Figure 2 and through an analogue 
VCO in Figure 1. The key function that allows tracking 
to take place within the data processor itself is the 
sliding discriminator. To form a frequency discriminator 
the power output of two neighbouring bins can be 
subtracted. If it is required to vary the position of 
a discriminator continuously then a technique of sliding 
the filters of the FFT filter will be needed to accomplish 
this. The following scheme is adopted. Given the 
complex output of two adjacent FFT bins it is possible 
to construct by using these two quantities alone a new 
FFT bin which will have as its peak any value of frequency 
between the mid-point of the two original bins. In order 
to derive a simple algorithm a rectangular window is 
assumed and phase factors between adjacent bins of T/N 
are ignored. The new FFT bin can be written approximately 
as 

f( x ) = * g ig * * ; p :sin^(l-x ) (2) 

where x is the fraction of a bin from the centre of the 
first FFT bin where the response is required. 
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lt can be shown that if a, 0 take the values 
- 6 2 D. + ucotit6 (1-6)3 



a = 



(3) 



Ticotn6 (6 (1-6) (1-26)) - 1 + 26 - 26 2 
P = 1 - a 

then the point f(x = 6) is a maximum of the filter of 
Equation (2). Hence a sliding filter can be implemented 
by taking the frequency f evaluating the bin nearest to 
this and then calculating the parameters a, 0 needed to 
give two filters separated by one half a bin separation 
from the target frequency (£)• The same technique is 
used on the filters regardless of the type of window 
that is being used. The equations for the sliding 
filters thus become :- 

N ta = Integral part of ( (f + 0.5B) t g + 1.5) 

x '» -0.5 + N. - N - f t 

ta p s 

= a(x) (where a is given by a pre calculated 

set of values) (4-) 

C a = C l b < N ta " 1) " (1 " V b(N ta ) 

C b = (1 - b(N ta + 1) - Cl b(N ta ) 

where is the doppler frequency according to the digital 
tracking loop, B is the PPT bandwidth, t is the update 
time, N p is half the number of bins in the PPT and b(N) 
is the N PPT bin (complex). Thus c and a are the two 
sliding filter outputs, and are a^ain complex. 



^SECRET 



\A< S-ECK &P 



1? 



-12- 

It is then necessary to construct the 
discriminator, and this can be done by taking 
d = | c al ~ | C b| • However the sliding filter technique 
results in a scaling at the origin of this quantity 
which depends on filter position and hence a correction 
for this has to be applied. A quadratic correction term' 
is used of the form 

S d = a x + a 2 x + a e 5E 2 (5) 

where x = min (x, 1-x) and a 1$ a 2> a^ take various values 
for the different window functions being used. This also 
allows the discriminator to be automatically scaled for 
the appropriate window function. 

The discriminator 'd* must also be automatically 
gain controlled and this is accomplished by a simple 
first order feedback system operating off a square law 
detector. The detector operates off a new bin of the form 

°agc - C a + °b < 6 > 
which gives a filter centred on the target frequency. 
In conjunction with the discriminator this is believed 
to give optimum performance.. The age is thus 

Agc Q = Agc n-1 k + (1-k) l c agc| 2 S c (7) 

where S c is another scale correcting factor for the 
different gains found at different positions of the 
sliding filter and different windows. This is given by 
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S c = 1 + 5c b 1 + x 2 b 2 (8) 
where b.^ and b 2 are constants depending on the type of 
window, k equals exp (-t H /t agc ) where t agc is the age 
time constant. The resultant discriminator output £ d is 

& d - d S c S d/ (A S c n V Hz (9) 
The digital tracking loop is closed entirely 
within the data processor (see Figure 2) and provides 
the estimate of target frequency f that is needed for 
the confirmation process. The particular implementation 
used is a type two loop. A block diagram of the system 
is shown in Figure 2. 

It can be seen that there are two outputs of 
the loop, one for the confirmation process (f) and one 
feedback (?) to the sliding discriminator. The transfer 
functions of these two are different and are given by 

£ _ z 2 ^ (1 + G) - gz 

£ + 2 (g (1 + G) -2) + (1 - g) 

~ 2 (10) 

| = + g (g — n 2 

z + z (g (1 + G) -2) + (1 - g) 

Where f is the target frequency relative to the centre 

of the IF sum channel. These transfer functions are 

implemented by taking the output of the discriminator £ * 

d 

and then multiplying by the gain g 1 (Gl). 



(M SECRET' 
-14- 

The algorithms for the implementation of these 
two closed loop transfer functions are as follows. If 
the agc'd discriminator output is £^ then 

f x (n) = + g x d 

f 2 (n) = f 2 (n-l) + g^Cn) 

to sliding f (n) = f,(n) + f p (n) (11) 

discriminator 

to confirmation ? (n) = f^n) + f 2 (n-l) 

The doppler loop described so far thus gives 
a signal f indicative of the target frequency f relative 
to the centre frequency of the speed gate (see Pig. 1), 
The frequency error £ d is added to f and multiplied by 
a variable gain G3 which feeds into the digital/analogue 
converter (DAC) of Figure 1. This signal labelled 
'DTL 0/P' is fed to an integrator to control a voltage 
controlled oscillator. Added in to the integrator input 
is a further input from an accelerometer to give a rate 
assistance component to the doppler tracking loop. 
The VCO signal is mixed with the LO signal and the. 
difference frequency is extracted in the single sideband 
mixer, amplified and transmitted from the aerial via the 
circulator. The loop settles when • DTL 0/P' is zero, 
i.e. when the transmitter frequency is such that the 
target return frequency is centred in the first speedgate(24-), 
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The parameters in the doppler loop are made 
adaptive to allow for changing target conditions so that 
optimum performance is maintained at all times. 
Referring to Figure 2: 

(i) Immediately after target detection, gains 
Gl, G2 and G3 are set to high values so 
that the loop settles quickly. The gains 
are then reduced until, after a pre- 
determined length of time, the gains take 
constant lower values to reduce the noise 
in the loop. The gains are increased as 
the missile flight progresses according 
to the signal strength measured in the age 
amplifiers (connection from age not shown), 
or by other signal to noise assessment 
methods (see later), and according to the 
FFT bin width selected (see below), 
(ii) During the missile boost phase, the doppler 
tracking loop experiences substantial 
transients and the FFT bin width is selected 
to be a large value (note the bin width 
equals the total FFT bandwidth divided by 
the number of FFT bins.) Gains Gl, G2 and 
G3 are set to high values to give a rapid 
response. 
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(ill) If a target at long range is to be acquired, 
the FFT bin width is set to be very narrow, 
and Gl G2 and G3 are set to low values after 
loop pull in. This optimises the acquisition 
threshold. 

(iv) If the target is at short range the PFT 
binwidth is set to be wide and Gl G2 and G3 
are set to high values. . This optimises 
acquisition speed and minimises tracking 
errors. A wide PPT binwidth gives a 
faster update rate and hence gives the 
optimum performance very dose to impact 
when rapid response is important. Target 
(2) range is determined from the range 
tracking loop (Figure 1). 

(v) Once the loop has pulled-in ,the signal f 
is proportional to the target frequency 
rate, i.e. the target acceleration. If 
f is large the PPT binwidth is set to a 
high value to allow a rapid tracking 
ability. The loss in doppler discrimination 
against other targets at similar speeds is 
not important since a manoeuvring target 

is not able to keep in close proximity to 
other targets. When ^ is small, a small 



,, r -\ - 




SECRET' 




FFT binwidth is used to give maximum 
doppler discrimination against targets 
flying in formation, 
(vi) When multiple targets in close proximity 
are detected (see later) or the target 
acceleration is small the FFT binwidth 
is set to a low value to give maximum 
doppler discrimination* 

The design of the overall angle tracking system 
can be conveniently divided into three parts. Firstly 
the EAT receiver, secondly the angle tracking filter 
and thirdly the servo control and stabilisation system. 

A block diagram- of the complete angle tracking 
and servo control system is shown in Figure 5. The 
angle tracking filter estimate of angle in the two 
channels (see below) is used to add a controlled portion 
of the sum channel which corresponds to these angles 
according to the stored D/S slope at the origin, into 
the difference channels as in Figure 1. These difference 
channels then pass through the same IF chain as the sum 
channel, are sampled and analogue -to-digital converted. 
2?hey are then Fourier Transformed to give the spectrum 
of signals on the azimuth and elevation channels. The 
sum and difference channels are agc'd digitally and 
elevation and azimuth phase-sensitive-detector outputs 
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are formed by taking the real part of DXS; these 
outputs are proportional to the angular error between 
the true target position and the best estimate of the 
target position. These calculations are all performed 
at the same update rate as the doppler and range tracking 
loops. 

After the PSD's it is necessary to perform an 
axis transformation in order to transform from receiver 
axes to line of sight axes. This is implemented digitally. 
These signals are then passed through an angle tracking 
(Kalman) filter which produces best estimates of sight 
line rate using a variable gain algorithm. The sight 
line rate signals are used in various ways; firstly they 
provide an output to the autopilot as the guidance 
command for which purpose they are transformed into 
missile axes; they are also used as the input to the EAT 
integrator which feeds estimate of angle back to the EAT 
element and finally it is used as input to the line of 
sight observer to drive the antenna mechanism. This 
latter is added to an inertial reference unit (IRU) gyro 
output in order to provide the total movement of the 
target as measured in line of sight axes. 

These angular rates are still calculated at 
the Q-Ow) update rate of the loop. These rates are then 
integrated using a high update digital integrator and 



M .SECRET 



^ - 



-19- 

and the output is transformed from line of sight axes 
into motor shaft angles using digital axes transformations. 
The angle demands are then used to control a position 
control loop for the antenna. The angle demand is 
subtracted from the potentiometer pick-off s and the 
error is fed into a controller which feeds demands to 
the antenna mechanism, servo and motors. The antenna/ 
reflector plate is thus position slaved to the integrated 
IRU gyro and the target sight line error. A high update 
rate controller is used. The effect of this is to remove 
body motions on the receiver output. An additional line 
of sight feed back is used to the EAT element. This is 
added to the other feedback path and ensures that the 
best estimate of sight line error can be used to subtract 
from the actual sight line error for EAT. The combined 
EAT feedback path has to be angle transformed back into 
receiver axes and a radome aberration correction is added 
in, as a function of the gimbal angles, as indicated in 
Figure 3. 

The angle tracking loop of Figure 3 will now 
be described in some detail. The sum and difference FFT 
outputs in Figure 1 are used to derive the target 
boresight error signal (i.e. the target angle off aerial 
boresight) in a manner illustrated in Figure 3. The 
sliding filters are similar to those described for 
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Figure 2 and are controlled to be centred at frequency f . 
The circuit forms the product of sum and complex conjugate 
of difference and takes the real part. This function is 
equivalent to that of a phase sensitive detector (also 
called a coherent detector). This signal is divided by 
the age detector output to give the amplitude ratio for 

each of the two channels i.e. (D_ - k S)S and 

e e 

(D — k S)S — *j 

a a where S , the mean square of the sum signal, 

is the output of the age detector. Only the elevation 
difference channel is shown but the azimuth difference 
channel is treated similarly* This is fed through a 
filter (e.g. a Kalman filter) to give an output £ . 

t is added to an aberration correction factor and 
fed to the EAT gain element to close the loop. The 
servo control output from the filter is a separately 
filtered version of 6 and drives the servo motors of 
the antenna mechanism. 

The aberration corrector gives outputs equal 
to the azimuth and elevation aberration levels of the 
radome as a function of the inputs. A function is 
stored in the memory of the corrector which takes 
account of the transmitter frequency selected, the 
measured temperature of the radome and the two gimbal angles, 
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When the EAT loop is in equilibrium £ gives the true 
measure of the boresight error, witty the aberration 
error subtracted . 



function of conditions during missile flight and in 
general are arranged to give high bandwidth when the 
estimated time to impact is small • Time to impact is 
estimated by dividing the missile-target range (from the 
range tracking loop) by the closing velocity (from the 
doppler tracking loop). The filter bandwidth is also 
made a function of signal to noise ratio, the bandwidth 
being low when the signal to noise ratio is low. 



ways. Firstly, the confirmation circuit of Figure 1 can 
be adapted to measure the S/N by forming the ratio of 
the power in the target FFT bin to the average power in 
the adjacent bins. Secondly, as shown in Figure 3, a 
signal to noise indicator can be derived by taking the 
imaginary part of. the product of the sum and difference 
channels and dividing by the age detector output. This 
gives an indication of the power of the signal components 
that are either incoherent between the sum and difference 
channels or in phase quadrature. The •incoherent. 1 power 
component is large when thermal noise or jamming in 
either the sum or the difference channels is large, and 



The Kalman filter parameters are varied as a 



Signal to noise ratio can be derived in two 
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tends to zero at high signal to noise ratio. The output 
is also large on extended targets or on targets flying 
in formation that are unresolved by the seeker. The 
ratio of the 'incoherent 1 power component to the signal 
to noise ratio indicated by confirmation-circuit is 
therefore an indicator of targets in close proximity, a 
large output indicating the presence of multiple targets, 
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